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There is an increasing interest in quantifying the predictive power in nuclear structure calcula- 
tions. We discuss how both experimental and systematic errors at the NN-level can be used to 
estimate the theoretical uncertainties by rather simple means and without solving the full nuclear 
many body problem. We emphasize the role of effective interactions defined by coarse graining 
the NN potential to length scales of the order of the minimal de Broglie wavelength probed be- 
tween nucleons in nuclei. We find an a priori error of AB/A ~ 0.1 — 0.4MeV for the binding 
energy per particle throughout the periodic table for 2 < A < 208, and a linear growth of the er- 
ror with density for nuclear matter AB/A ~ 3.75p n . m . an d neutron matter AB/N ~ 3.5p n .. This 
suggests to limit the computational effort in solving the Nuclear Many Body Problem to such an 
accuracy. 
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1. Introduction 

The present contribution is based on our recent work [|TJ, where an old issue is adressed from 
a modern perspective. The key question is, can we quantify errors in binding energies based on our 
incomplete knowledge of the NN interaction ?. 

The nuclear many body problem consists of diagonalizing the Hamiltonian H with multinu- 
cleon interactions, 

H =L T i+L v 2,ij+ L v xijt+ £ v Am +... (i.i) 

i i<j i<j<k i<j<k<I 

where the indices run 2 < i < A and 7} = pf / (2M N ) is the kinetic energy of the i— th nucleon, V%jj 
the NN potential, V^jjt the NNN potential, etc. Roughly speaking one can fix V2 from the deuteron 
and NN-scattering data, V3 from triton and nucleon-deuteron scattering, V4 from the a— particle 
and any scattering process involving four nucleons, etc. 

First principles calculations in Nuclear Physics have traditionally been dominated by the idea 
that once the fundamental NN-interaction is accurately known one is left with the intricacies and 
complexities of the many body problem on the theoretical side. Actually, ab initio calculations 
in Nuclear Physics have been carried out as if there was complete knowledge on the elementary 
NN-dynamics. Therefore, a long-term effort has been carried out to substantiate this assumption 
by continuously improving the NN-potentials, but little attention has been paid to determine the 
uncertainty in the potentials themselves. We have recently filled this gap by carefully analyzing 
different error sources ^ and using the concept of a coarse grained potential [g]. 

2. NN potentials 

Along the years many studies have been oriented towards improving the NN-potentials by 
performing # 2 -fit of a Partial Wave Analysis to the abundant available pp and np experimental 
data. From a purely statistical analysis point of view, we remind that in order to determine reliable 
confidence levels on the fitting parameters, i.e. error estimates, one must have # 2 /d.o.f < 1. Larger 
values of the reduced % 2 (most Bonn,Paris, Esc or Nijm93 potentials produced # 2 /d.o.f ~ 2) would 
actually diminish the uncertainties, as there is a large penalty to change the most likely fit parameter 
values, and would produce unreliable error estimates. Thus, an error analysis of NN phase-shifts 
for several partial waves became first possible when the Nijmegen group [Q] carried out a Partial 
Wave Analysis (PWA) fitting about 4000 experimental np and pp data (after rejecting further 1000 
of 3a-mutually inconsistent data) with # 2 /dof ~ 1. 

In general, fits consider differential cross sections, polarization asymmetries, etc. for a given 
set of energies and angles (£,-, 0,). While the partial wave expansion allows to evaluate any scat- 
tering angle, the fixed energies allow to determine phase shitfs only at those measured energies, so 
that the phase-shifts themselves become independent fitting parameters. On the other hand, the an- 
alyticity of the S -matrix based on the meson-exchange picture, guarantees that the phase-shitfs are 
smooth functions for real positive energy. A handy way to combine this energy-angle information 
and implementing the expected smooth energy dependence is by introducing an auxiliary potential 
with a set of adjustable parameters. The bench-marking Nijmegen fit fixed the form of the poten- 
tial to incorporate charge dependence. It contains an energy dependent square well operating at a 
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distance below 1.4fm, a One-Pion-Exchange (OPE) contribution starting at 1.4fm, a One-Boson- 
Exchange (OBE) piece below 2 — 2.5fm and an electromagnetic contribution. Energy dependence 
reflects retardation in the interaction due to unobserved virtual excitations. However, it is incon- 
venient to perform Nuclear structure calculations since it requires solving a time-dependent many 
body problem. At present there are a variety of NN (energy independent) potentials fitting a large 
body of scattering data with # 2 /dof ~ 1 [Q, ||, ^, ^, ||] which allow the application of conventional 
stationary dynamics. Surprissingly, error estimates on the potential fitting parameters are never 
provided. In [Q] we analyze different error sources, in particular the short distance non-localities 
of the interaction which may depend on energy, linear and angular momentum. While in principle 
these non-localities are on-shell equivalent (see e.g. Ref . [Q] for a proof in a 1 /Mn expansion) they 
generate quantitative differences. 
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Figure 1: np scattering observables for different energies in the laboratory system as a function of the CM 
angle. Left panel: Polarization at £lab = 325MeV. Middle panel: Differential cross section at £lab = 
324.1MeV. Right panel: D t depolarization at £lab = 325MeV. The short-dashed line denotes our results 
by the coarse grained interaction [^]. The band represents the compilation of the PWA and five high quality 
potentials [0,^,0,^] which provided a X 1 Id.o.f < 1. For references for the experimental data see http://nn- 
online.org and http://gwdac.phys.gwu.edu/. 

From our reanalysis of the discrepancies in the different partial waves and have found that 
the PWA statistical errors quoted in the original Nijmegen work [ffl|, A5p WA , turn out to be smaller 
than the discrepancies among the fitted phases stemming from the different potential fits, separately 
claiming a # 2 /d.o.f. < 1 to data. In other words, 



A<5pWA ^ y^3jI>»- 5 ) 2 > $i — ^tijmn,^leid93)^rijmI)5Avl8)5speC)5cDBomi (2-1) 

Thus at £lab = 350MeV in the 'So channel one has a s.d. of 0.7° whereas A5pwa = 0.3°. This 
counter-intuitive result relies not only on the specific forms of potentials which treat the mid- and 
short-range behaviour of the interaction differently but also on the fact that the fits are mainly done 
to scattering amplitudes rather than to the phase-shifts themselves. The systematic discrepancies 
are vividly illustrated by looking at Fig. [j] where the spread of the different potentials and the PWA 
analysis increases in the region where no data constrain the result. 
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3. Direct error estimates 

The most direct way of quantifying binding energy uncertainties would be to undertake large 
scale ab initio calculations using the different two-body potentials, say V 2 with i = l,...N, yield- 
ing flW (A) whence a mean B(A) and a standard deviation Afi(A) can be constructed, 



5 ( A ) = AB(A) = 0J_£(HW(A) -B(A))2 (3.1) 

For instance, the triton binding energy obtained by Faddeev calculations is 8.00, 7.62, 7.63, 7.62, 7.72 



and 8.50 MeV for the CD Bonn [jlOJ, Nijm-II, Reid93, Nijm-I, AVI 8 [jll|] and the covariant spec- 
tator model ^ respectively. This yields the combined result B^ = 7.85(34)MeV (exp. B3 = 
8.4820(1 )MeV) i.e. Afi 3 /3 = O.llMeV. Note that a three-body interaction would account for 
the missing lMeV-binding. 

This error estimate procedure stops beyond the A = 4 nucleus, due to computational and the- 
oretical difficulties related to the form of the potential. From an ab initio viewpoint, only Monte 
Carlo calculations may go up to A = 10 when potentials are fixed to be r— dependent with a nonlo- 
cality in terms of the relative angular momentum operator, becoming the standard Hamiltonians for 



ab initio calculations of light nuclei [ |12| ] and dense matter [|13|]. Again, a lack of error estimate of 
the potential parameters makes it impossible to deduce the theoretical errors in such calculations. 
On top of that, one should add the systematic errors due to the different forms of the potential. 

Finite experimental accuracy provides also an statistical error AV n on the n-body force. In 
addition, unless the most general possible n-body forces are considered some systematic error is 
introduced. On the other hand, the definition of the three- and higher-body interaction depends 
on the two body potential, so any uncertainty in the two-body interaction will carry over to the 
three-body interaction. Thus, even if we fix it say in the A = 3 system, there will always be a 
residual uncertainty in the A + 1 = 4 calculation. To fix ideas it is convenient to think of the n- 
body potentials as random variables. They depend on some unknown parameters (ci , . . . c/v) and 
will eventually be determined from # 2 -fit to some data involving n-body interactions. This would 
provide a probability distribution of parameters P(ci,...cn), according to which the normalized 
variable v n = (V„ — (V,,))/ \J ((V n — {V n )) 2 ) which has zero mean, (v n ) = 0, and unit variance, Av„ = 
1, can be defined. Correlations between the two- and three-body forces can be ignored within errors 
if ((V2V3)) 2 < (((V2V3 — (V2V3}) 2 }. Although establishing the validity of this condition requires a 
thourough analysis of multinucleon forces, we will assume this to be the case, so that 

AB(A) 2 = AV 2 2 + AV 3 2 + AV 4 2 + . . . (3.2) 

in which case, estimating the two-body uncertainty provides a lower bound on the total uncertainty. 
Note that non-vanishing statistical correlations may change this. 

4. Error estimates from Coarse grained interactions 

In our recent work [Q] we have shown how coarse grained interactions [^] can advantageously 
be used for the purpose of error estimate. Besides making accurate fits to the combined NN- 
scattering database obtained from the PWA and the 6 high-quality potentials this representation 
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allows to side-step the short distance complications arising in nuclear structure calculations. This 
is very much in spirit of the Vi ow k-approach [14], the Unitary Correlation Method (UCOM) [|5|] or 
the Similarity Renormalization Group (SRG) (see e.g. jl6| ] and references therein), where smooth 
soft-core potentials are produced after the original potential has been evolved to the relevant scale. 
However, these transformations are applied to the existing high-quality potentials and become com- 
putationally cumbersome. We suggest instead to determine the effective coarse grained interactions 
directly from fits to the scattering data up to a given maximum energy. For light nuclei A < 40 
it turns out that a reasonable maximum CM momentum is about pcm ~ 200MeV ^ since the 
highest wavelength resolution l/p ~ lfm is much larger than the short range repulsion distance, 
flcore ~ 0.5fm, present in central waves. Thus, we expect mean field calculations to be applicable. 

In this way the systematic uncertainties due to the different forms of the potential can be 
propagated. We have found that for doubled-closed shell nuclei S-waves dominate the uncertainty. 
Using Harmonic Oscillator Shell model wave functions we get find [p 

=0.08(1), ^£=0.12(1), =0.28(2), = 0.34(2)MeV (4.1) 

3 4 16 40 

the errors depending on the fitting cut-off LAB energy. While one may reasonably doubt that 

variations in the binding energies can be monitored by coarse grained interactions with harmonic 

oscillator wave functions, we have found [P that this is not the case by analyzing variations of 

binding energies in ab initio calculations induced by a relative 1% change in the AVI 8 potential 

parameters [|l7|]. 

For heavier nuclei a simple estimate of errors due to the two body interaction uncertainty 
can be made by using Skyrme effective interactions (for a review see [18]) where for symmetric 
N = Z = A/2 nuclei one has 



AB 3 , 

— = — M) 
A 8A 



Jd 3 xp(x) 2 , (4.2) 



where to = J d 3 x(Vi So + Vi Sl )/2 and xoto = J d 3 x{— Vi So +Vi Si )/2. Propagating errors we get to = 
— 0.92(1 )GeVfm 3 , in agreement with the equation of state used by the Trento group @ at low 
densities, to ~ — 0.9(1 )GeVfm 3 and with coarse graining of NN interactions in CM momentum 
space down to A ~ 0.3GeV gives a compatible value, to ~ — 4tt 2 / (M#A) p0[]. We may implement 
finite size effects by using a Fermi-type shape for the matter density p(r) = po/(l + e( r ~ R )/ a ) 
with R = roA? and r = l.lfm and a = 0.7fm and normalized to the total number of particles 
A = J d 3 xp(x) we get a result compatible with Eq. (1.1) and depending on the value of A for 
4 < A < 208. In Fig. @ we illustrate the situation by imposing our error band on the binding 
energy of stable nuclei. For nuclear and neutron matter the effect grows linearly with the density 
(in fm~ 3 ) 

= ht p ~ 3.75p , ^ = U[t (\-x )]p n ~3.5p n (4.3) 

In Fig. |2| we implement our error estimates on the EOS calculations of the Trento group [|^], where 
their uncertainties reflect the accuracy in solving the many body problem only. As we see, their 
errors are much smaller than those estimated here. We remind that among other effects we have 
neglected possible statistical correlations between 3A^ and 2Af forces, thus a more thorough analysis 
incorporating n-body coarse grained interactions would provide a more definite answer. 
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Figure 2: Superposition of estimated errors based on propagating Nucleon-Nucleon uncertainties via coarse 
grained interactions. Left panel: Stable nuclei. Middle panel: Nuclear Matter Equation of State (EOS). Right 
Panel: Neutron Matter Equation of State. EOS Data are from the Green Function MonteCarlo Calculation 
of the Trento group [|l9j, 
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